Carrier localization mechanisms in InGaN/GaN quantum wells 
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Localization lengths of the electrons and holes in InGaN/GaN quantum wells have been calcu- 
lated using numerical solutions of the effective mass Schrodinger equation. We have treated the 
distribution of indium atoms as random and found that the resultant fluctuations in alloy concen- 
tration can localize the carriers. By using a locally varying indium concentration function we have 
calculated the contribution to the potential energy of the carriers from band gap fluctuations, the 
deformation potential and the spontaneous and piezoelectric fields. We have considered the effect 
of well width fluctuations and found that these contribute to electron localization, but not to hole 
localization. We also simulate low temperature photoluminescence spectra and flnd good agreement 
with experiment. 



Despite the relatively high dislocation density in In- 
GaN/GaN quantum well (QW) structures grown on c- 
plane sapphire, the room temperature photolumines- 
cence (PL) efficiency of such QWs can be very higtP. 
This is generally attributed to localization of the carriers 
which reduces the effect of non-radiative recombination 
at the dislocationPHS, xhe precise nature of this local- 
ization is still however a matter of debate. Three possi- 
ble causes of carrier localization that have been widely 
cited are well width fluctuationJ^HZl^ random alloy fluc- 
tu ationa ^ including In-N-In chain^, and indium cluster- 
in ^io | ii | Smeeton et al. have shown however that gross 
indium clustering observed in TEM images of InGaN 
QWs can be caused by electron beam damage^^. 

In this work we present the results of theoretical cal- 
culations which demonstrate the importance of random 
fluctuations in alloy composition, and of well width fluc- 
tuations. Previous theoretical work on carrier localiza- 
tion has considered the properties of bulk InGaN with 
or without embedded InGaN quantum dots (QDsJP^Ell. 
Using an atomistic empirical pseudo-potential method 
they found that in bulk zinc-blende InGaN hole wave- 
functions are localised by alloy fluctuations alone, and 
that the electron wave-functions required a dot, or clus- 
ter like confinement to be localised. The properties of 
these bulk zinc-blende materials are expected to differ 
from those in the wurtzite QW structures we consider 
because of the effect of the large strain, and piezoelectric 
field across the QW in the case of wurtzite materials. 

Our method is to calculate the potential energy land- 
scape of the QW using a random distribution of indium 
atoms. Once we have calculated this potential energy we 
solve the effective mass Schrodinger equation using the 
finite difference approximation to find the energy eigen- 
values and wave- functions of the carriers. Performing this 
calculation for different distributions of indium atoms 
and averaging the results reveals the mean localization 
lengths for electrons and holes. 

Experimental results showing the random nature of the 



indium distribution and the occurrence of well width fluc- 
tuations are described in Section |ll Our calculations of 
the carrier potential energy and the subsequent solution 
of the Schrodinger equation are described in Section [llj 
and the results of these calculations are discussed in Sec- 
tion |TlIl 



I. MICROSTRUCTURAL DATA 

In order to model the localization of carriers in InGaN 
QWs, a realistic description of their nanoscale structure 
is required. Here, we base our model on the experimen- 
tal data of InGaN QWs gained by Galtrey et al.E^ us- 
ing atom probe tomography (APT). APT has two ma- 
jor advantages over conventional transmission electron 
microscopy (TEM) studies of InGaN QW microstruc- 
tures: firstly, the three-dimensional nature of the APT 
data set is more appropriate to the development of a 
three-dimensional model than the two-dimensional pro- 
jections usually recorded in TEM, and secondly the high 
energy electron beam used in TEM has been reportecP^l 
to rapidly damage InGaN QWs, resulting in potentially 
unreliable structural data. We will outline the key re- 
sults of APT studies of InGaN quantum wells grown at 
a single temperature, and highlight the specific features 
of the APT dataset that we have used as an input to our 
model. 

For quantum wells grown at a single temperature, with 
a; = 0.18 and 0.25, the APT data revealed the In^^Gai-^N 
within the quantum well to be a random alloj^^^j xhe 
experimentally-determined distribution of indium atoms 
within the quantum well was compared to the expected 
binomial distribution for a random alloy. Figure [l] shows 
this comparison for a typical x — 0.25 quantum well sam- 
ple. A test was used to assess whether there was any 
statistically significant deviation from the random distri- 
bution, but no evidence was found for any indium cluster- 
ing^. Hence, in our calculations, the InGaN within the 
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FIG. 1: Statistical analysis of the distribution of indium con- 
tents in InGaN QWs indicates no deviation from that ex- 
pected in a random alloy. For the composition frequency dis- 
tributions shown here, APT data from green-emitting QW 
material was divided into bins of 100 atoms, and the indium 
content was calculated for each bin. The solid line shows the 
binomial distribution that would be expected in the case of a 
random alloy. Similar data are also available in Ref^. 



QW has been modeled as having a random, uncorrelated 
indium distribution. 

Further APT studieJ^^ revealed important information 
regarding the two QW interfaces. We describe the in- 
terface in which the InGaN is grown on the GaN as 
the lower interface and the interface at which the In- 
GaN is capped with more GaN as the upper interface. 
The lower interface was found to be both smooth and 
abrupt. In contrast, the upper interface was both more 
diffuse and rougher^^. In indium iso-concentration sur- 
faces illustrating the upper interface, such as that shown 
islands one to two monolayers high and a 



in Fig 2 (a 



few nanometers across were observed on the InGaN sur- 
face, and this observation was supported by other data 
from TElvP and atomic force microscopy (AFM)^^. In 
our model, we have included islands at the upper inter- 
face similar to those observed experimentally, since the 
well width fluctuations (WWFs) the islands cause have 
previously been predicted to significantly affect carrier 
localizatiorPt^. We have investigated the effect of lOnm 
diameter, one mono-layer thick disk-shaped WWFs in 
our c alculations, as schematically demonstrated in Fig. 
2(b) We have not included interface diffuseness in our 



model, but since this diffuseness is homogeneous in the 
plane of the QW we expect this to have a less pronounced 
effect on the localization than the WWFs. 




FIG. 2: (a) An iso-concentration surface based on APT data 
illustrating the roughness of the upper interface of an In- 
GaN QW. Islands one or two monolayers high and a few 
nanometers across are observed-'^, (b) Schematic diagram 
of InGaN/GaN QW with disk-shaped monolayer well width 
fluctuation. 



II. THEORY 

To investigate the localization of the carriers in these 
structures the carrier wave functions are calculated by 
solution of the effective mass Schrodinger equation. This 
method requires knowledge of the potential energy for 
the carriers. Our calculation of the potential energy of 
the carriers includes the effects of the band offset between 
InN and GaN, the spontaneous polarization, the piezo- 
electric field, and the deformation potential. As this po- 
tential varies spatially in all three dimensions, this poten- 
tial energy can be referred to as an "energy landscape" . 
We discuss the calculation of the potential energy in sec- 
II A| and the solution of the Schrodinger equation 
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A. Potential Energy Landscape 

The potential energy landscape depends critically on 
the distribution of indium atoms. As described above 
we take the indium atoms to be distributed randomly in 
the QW. For the continuum strain mechanics we employ, 
the indium distribution must be described by a smoothly 
varying concentration function. To calculate this concen- 
tration function we first occupy cation lattice sites with 
indium atoms with probability equal to the nominal con- 
centration. We use Gaussian functions with a standard 
deviation equal to the cation lattice spacing to smooth 
this positional data. An example of the resulting indium 
concentration function, x(r), is shown in Fig.js] The fiuc- 
tuations in indium fraction seen in Fig. [3] arise naturally 
from the random distribution of indium atoms and is fun- 
damental to the calculations which follow. It should be 
noted that the maxima of the indium concentration seen 
in Fig. [3] are not due to individual indium atoms, but 
instead represent the averaged effect of several indium 
atoms lying close to one another; this can be confirmed 
by examining the distance scales on the axes of the plot. 
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FIG. 3: A result of our calculation of the indium fraction in 
a QW of width 3nm with an average indium concentration of 
25%. The z direction is normal to the x-y plane of the QW. 



These fluctuations in local indium fraction are respon- 
sible for the local variations in quantum confinement. 
To calculate this locally varying band gap we use the 
quadratic approximation, 



GaN 



(1) 

Eg"^ , and E^°-^ , are the energy gaps in InN and GaN re- 
spectively, and b is the bowing parameter. These parame- 
ters and all the relevant mechanical constants are taken, 
or derived, from Vurgaftman and Meyer's review arti- 
clj^. and the effective mass parameters are from Rinke 
at al.Csl. 

The strength of the spontaneous polarization which 
occurs in wurtzite structures also depends on the indium 
concentration and we describe the polarization using a 
linear approximation. 



P,,(r)=x(r)Pi;'^ 
We can also write this as 
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As the polarization is purely in the z direction, we can 
write Pl'^^ - P^jf ^ = APz. This can be used in Pois- 
son's equation. 



Psp = -AP-^ = -V • (eoerV0,p). 
oz 

which, after taking the Fourier transform, gives 



(t^spCk.) 
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for the corresponding electrostatic potential. This contri- 
bution is evaluated numerically and transformed to real 
space using a fast Fourier transform (FFT). 



Further to this the lattice constants of InN and GaN 
differ, so that the fluctuating indium distribution causes 
an inhomogeneous strain field. This strain causes a defor- 
mation potential and a large piezoelectric field. We use a 
Green's function method to find the strain field, and the 
following description closely follows that by Andreev et 
al.SH. The Green's tensor Gfc„(r) gives the displacement 
at r in the direction fc, due to a unit point force in direc- 
tion n placed at the origin. For an infinite anisotropic, 
homogeneous elastic medium, it is the solution to the 
equatiorPSI 
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For small deformations we are then able to write an ex- 
pression for the strain tensor in terms of the Green's func- 
tion and concentration function, x{^)- We solve for the 
Green's function by Fourier transform and obtain an ex- 
pression for the strain in Fourier space. By numerical 
evaluation of this function, and by performing a FFT we 
are able to calculate the real space strain distribution due 
to the inhomogeneous indium concentration function. 

The deformation potential accounts for the change in 
energy levels due to this strain. To calculate the defor- 
mation potential the band gap energy is given in terms 
of the deformation potentials, D and a, and the strain 
tensor U by: 



Eg{\j) = Eg{Q) + U,, (ai - D-i) + {U^^ + Uyy) (fls 



^ (7) 
where the r dependance is implicilpD. 

Once we have calculated the strain in the material find- 
ing the piezoelectric field is straightforward. The polar- 
ization caused by the piezoelectric effect can be written 
as Pi = Ci^kiUki, where Ci^ki is the piezoelectric tensor. 
Poisson's equation is solved in the same manner as for 
the spontaneous polarization. 



B. Carrier Wavefunctions 

The carrier wave functions and energy eigenvalues are 
obtained by solution of the effective mass Schrodinger 
equation. 



1 



m(r) 



V^(r) +V{r)^j{r) = E^P{r), 



(8) 

where we use an anisotropic, inhomogeneous effective 
mass. The hole in-plane (perpendicular to the c-axis) 
effe ctive mass increases rapidly away from the band 
edg^^mH gj^^ ^ Y) calculations for the c-plane QW sam- 
ples detailed below revealed holes ^ 15 — 55meV away 
from the band edge. To avoid the use of prohibitively 
expensive iterative methods we have assumed far from 
the band edge hole effective masses as given by Chuang 
and Chang2^. 
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III. RESULTS AND DISCUSSION 
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FIG. 4: The calculated ground state probability density of an 
electron (red) and a hole (green). The iso-surfaces are taken 
at one standard deviation of the maximum of the function, 
and plotted in real space in units of nanometers. Figure (a) 
shows the view along the z-axis, and (b) shows the view along 
the y-axis. The plane of the QW is normal to the z-axis, and 
the interfaces are a.t z = 3.5 and 6.5nm. 



As our potential energy V{r) has already been evalu- 
ated numerically, it is a straightforward exercise to solve 
this equation using a finite difference method, with pe- 
riodic boundary conditions on the wavefunction. The 
resulting sparse matrix eigenvalue problem is solved us- 
ing the ARPACK subroutine librarji^Sl. This allows us to 
find the lowest energy eigensolutions of the system. 

In Fig. |4]we show the calculated ground state probabil- 
ity density of an electron (red) and a hole (green); they 
were calculated for a 3nm wide 25%-indium QW, and 
are clearly localised on a scale of just a few nanometers. 
While Fig. |4] shows both hole and electron wavefunctions 
it should be noted that the interactions between them is 
not included in these calculations. To quantify this lo- 
calization we define our localization length as Ar, using 
the standard quantum mechanical expression 



(Ar) 



|(r-(r))^| 
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(9) 



where 



(10) 



To investigate the localization properties of the car- 
riers we performed the calculations described above for 
different random distributions of indium atoms. Calcula- 
tion system sizes of 15 x 15 x 10 nm were used to find the 
15 lowest energy hole states for 45 different distributions. 
The four lowest energy electron states were calculated in 
65 larger 45 x 45 x 10 nm systems owing to the larger 
localization lengths of these states. The systems stud- 
ied were 3nm thick QWs which varied in average indium 
composition between 5% and 25%. We investigate the 
effect of the local fluctuations in indium concentration in 
polar (c-plane) QWs, monolayer WWFs as described in 
section [ij and matching non-polar (a- or m-plane) sam- 
ples. 

The mean localization lengths Ar, as defined in Eq. 
|9] for the lowest energy hole and electron eigenstate s are 
plotted against the nominal indium fraction in Figs 5(a) 
and |5(b) respectively. Considering the effects of indium 
fluctuations only, the hole wavefunctions show a constant 
localization length of ~ l.Onm across the range of indium 
concentrations. The electrons show an increasing local- 
ization length from ~ 7nm for the 25% indium sample, to 
^ lOnm for the 5% indium sample. This trend of increas- 
ing localization length with decreasing nominal indium 
concentration is caused by the reduction in magnitude 
of the fluctuations (strictly the standard deviation of the 
concentration distribution decreases) as the nominal con- 
centration gets smaller. The holes do not display this 
behavior as the larger hole effective mass means they are 
less sensitive to the changing depth of fluctuations. The 
difference in overall localization length-scale can also be 
attributed to the difference in the effective mass between 
the electrons and holes. 

When we consider the added contribution of WWFs 
on the carrier localization lengths we see a marked dif- 
ference in their effect on electrons and holes. The mean 
localization length for the electrons in the 25% indium 
sample is reduced by a third with the introduction of 
this WWF, whereas for the holes the WWF has almost 
no effect on the localization. This difference arises from 
the effect of the electric field, which causes the holes to be 
localised on the smoother (lower) interface, and the elec- 
trons to be localised on the rougher (upper) interface. It 
should be noted that the holes are localised on a shorter 
length scale than the WWF. The effect of the large elec- 
tric field across the QW was investigated by comparing 
the polar results with calculations on QWs in a non-polar 
(a- or m-plane) orientation. For the polar QWs we have 
neglected the excitonic interaction, which is justified in 
wide wells due to the separation of the carriers caused by 
the electric field. This is not the case in non-polar QWs; 
nevertheless the calculations provide at least qualitative 
information about carrier localization in such structures. 
We also only consider the effect of alloy fluctuations in 
this orientation, not WWFs. We see that again the hole 
localization lengths are unchanged throughout the range 
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(a)Ground-state electron localization lengths 



10 



Without WWFs + 

With WWFs X 

Non-Poiar 

Experiment □ 



0.3 




0.05 0.1 0.15 0.2 0.25 

Average nominai indium fraction 



0.3 



FIG. 6; The mean indium fraction within one standard de- 
viation of the maximum of the electron and hole probabiUty 
densities is plotted against the nominal indium fraction of the 
QW. Error bars represent the standard error in the mean. 
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(b)Ground-state hole localization lengths 

FIG. 5: Mean ground-state localization lengths are shown for 
QWs of different indium fraction, with and without WWFs 
and in a non-polar orientation. Localization data from Gra- 
ham et al^ are also plotted for comparison. Error bars rep- 
resent the standard error. 



of concentrations investigated. The electron localization 
is also unchanged, suggesting the electric field has lit- 
tle effect on carrier localization. At this stage we have 
not yet considered the effect of the electric field on elec- 
tron localisation in the presence of WWFs. We have 
also compared our results with the localization lengths 
of Graham el. al^ in Fig. 5(a) They obtained Huang- 
Rhys factors from the strength of the LO-phonon replicas 
in PL spectra at low temperature and compared these 
with the results of calculations which assumed Gaussian 
wave functions for the localised carriers. This allowed 
them to estimate the carrier localization lengths. These 
show encouraging agreement with our hole localization 
results, which may demonstrate that the relevant local- 
ization length being probed in the experiment is that of 
the holes. 



shows the mean indium fraction in the neighborhood of 
the localised carriers for each nominal concentration. We 
are able to calculate this by comparing the local concen- 
tration function for the sample with the calculated wave 
functions, and take the neighborhood of the carrier to be 
within one standard deviation of the peak of the wave 
function. Figure |6] clearly indicates that the holes are 
localised in naturally occurring regions of higher indium 
concentration. In the nominally 25% sample for example, 
the mean indium concentration under the lowest energy 
hole wave functions was 37%. It should once again be 
noted that these above average indium regions occur due 
to random fluctuations in composition. The mean in- 
dium concentration around the electrons in this sample 
was approximately 19%. This is lower than the nominal 
concentration as the electron wave functions penetrate 
into the barrier. 

We can approximate the localization energies of the 
carriers by fitting the density of localised states with 
stretched exponentials. 



9e/h{E) « noexp 



\E-E, 



e/h\ 



(11) 



Considering only the polar QWs with no WWFs, Fig[6] 



where a = 1.2 for the holes and a = 2 for the electrons. 
We find the localization energy (cr) of the holes in the 25% 
indium c-plane sample without WWFs to be ~ 40meV, 
and ~ 13meV for the electrons. This again emphasizes 
the stronger localization of the holes in these structures. 

All of our results point towards a picture of localiza- 
tion in which the holes are strongly localised by statis- 
tical fluctuations in indium concentration, and electrons 
are localised on a longer length-scale by a combination 
of these alloy fluctuations, and well width fluctuations. 
However including the coulombic interaction may change 
this picture, as the strongly localised holes may "bind" 
electrons to them. 
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FIG. 7: A simulated PL spectrum for a 25% indium single 
QW with no WWFs, compared with low temperature exper- 
imental PL spectra reported by Graham et al.^. The low en- 
ergy feature of the experimental PL spectra are LO-phonon 
replicas of the main peak; the electron-phonon interaction is 
not included in the simulated spectrum. 

Our calculations allow us to determine the oscillator 
strengths and energies of photons emitted during recom- 
bination. To ensure that we included all of the localised 
hole states within a given sample, we performed the cal- 
culation in a slightly smaller system, and calculated more 
states. We calculated 4 of the lowest energy electron 
states in 50 48 x 48 x 10 nm samples, and 35 hole states 
in 800 12 X 12 x 10 nm samples. We approximate the low 
temperature PL spectra of the samples by assuming that 
the localised states are all occupied by carriers, and by 
convolving the discrete recombinations with Gaussians 
of standard deviation 7meV. This width is chosen to be 
small enough for it not to affect the width of the final 
spectrum, but large enough for the spectrum to be rel- 
atively smooth. Note that due to the finite size of the 
system there will be long wavelength fluctuations in the 
indium concentration which are neglected. The effect of 
these would be to broaden the low energy tail slightly. 
Comparison of our calculated PL spectrum for a single 



3.3nm thick 25% indium QW with no WWFs with an 
experimental PL spectrum of a nominally identical sam- 
ple is shown in Fig. [7j The width of the peaks, often 
attributed to localization effects, compare very well: the 
simulated PL has a FWHM of 69meV, and the FWHM 
of the experimental PL is 63meV. In our calculations 
the strong hole localization is the cause of this broad- 
ening. The low energy features in the experimental PL 
are LO phonon replicas of the main peak and do not ap- 
pear in the simulated spectrum as we have not included 
the electron-phonon interaction in our calculations. In- 
cluding the effects of WWFs increases this broadening on 
the low energy side of the peak, but the effect is much 
smaller than that caused by the hole localization. With 
allowance for the uncertainty in the experimental deter- 
mination of, for example, the piezoelectric coefficients, 
the peak energies also agree reasonably well. 

IV. SUMMARY 

We have been able to build a picture of the localiza- 
tion of the carriers in InGaN/GaN QWs by performing 
effective mass calculations for different random distribu- 
tions of indium atoms. We have found that the holes 
are strongly localised in regions of above average indium 
content. The electrons are less strongly localised by the 
indium fluctuations, but become more localised by mono- 
layer well width fluctuations. Our work demonstrates 
carrier localization without the need for gross indium 
clusters, which goes some way to explaining the carrier 
localization mechanisms which are so important for the 
operation of InGaN optoelectronic devices. 
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